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Ratios of azimuthal angle correlations between two jets produced at large rapidity separation 
are studied in the JV = 4 super Yang-Mills theory (MSYM). It is shown that these observables, 
which directly prove the SL(2, C) symmetry present in gauge theories in the Regge limit, exhibit 
an excellent perturbative convergence. They are compared to those calculated in QCD for differ- 
ent renormalization schemes concluding that the momentum-substraction (MOM) scheme with the 
Brodsky-Lepage-Mackenzie (BLM) scale-fixing procedure captures the bulk of the MSYM results. 

1 Introduction 

In recent years there have been many efforts to unveil the structure of the Af = 4 super Yang-Mills 
(MSYM) theory, sometimes considered as "the harmonic oscillator of the 21 st century" 1 since 
it seems to provide, in the planar limit, the first example of a solvable nontrivial quantum field 
theory in four dimensions. The discovery of an unexpected simplicity in the structure of scattering 
amplitudes [2-6] and the connection with string theory provided by the anti de Sitter / conformal 
field theory (AdS/CFT) duality [7-9] support this reasoning. 

In the present work MSYM is taken as a theoretical laboratory to compare with observables cal- 
culated in Quantum Chromodynamics (QCD). This is useful since, for scattering amplitudes at high 
energies, both theories are equivalent at the leading order (LO) level, with the higher loop results 
in MSYM being, due to the enhanced symmetries and particle content, much easier to compute. 
As a very interesting connection between both theories, it is known that the MSYM amplitudes, 
interpreted as part of QCD amplitudes and contributing to anomalous dimensions of gauge-invariant 
operators, provide the "highest degree of trascendentality" terms [10]. This property is now a key in- 
gredient to calculate anomalous dimensions to all-orders using integrability -based approaches [11]. 

'Preprint numbers: LPN1 1-30, IFT-UAM/CSIC-1 1-38, FTUAM-1 1-48 

'In a similar fashion, black holes can be thought of as the "hydrogen atom" of quantum gravity [1]. 
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Different observables have been proposed in the literature to be calculated in the context of 
MSYM [12-15]. The energy flow functions defined in terms of correlation functions of the energy- 
momentum tensor in [16] is a very interesting example. In QCD, typical observables are inclusive 
jet cross-sections. In the present work the study of ratios of azimuthal angle correlations in inclusive 
dijet production, for configurations where the two tagged jets are well separated in rapidity from 
each other, is carried out in detail. These ratios were first introduced in [17, 18]. 

In SECTION 2 the kinematical set-up under study is introduced and the calculation of the in- 
clusive dijet production in the Balitsky-Fadin-Kuraev-Lipatov (BFKL) framework [19-21] at next- 
to-leading order [22, 23] (NLO) both in QCD (already known [17, 18, 24]) and MSYM (a new cal- 
culation) in the minimal subtraction renormalization scheme (MS) is performed. In SECTION 3 
the momentum-subtraction (MOM) scheme [25-27] with the Brodsky-Lepage-Mackenzie (BLM) 
scale-fixing procedure [28] is applied to the QCD calculation with the intention to challenge the 
common statement that it captures the bulk of the conformal contributions to all orders in the cou- 
pling. If this is correct it should then give a very similar result to that previously obtained in MSYM. 
In SECTION 4 it is shown that this is indeed the case, in particular, for the ratios of azimuthal angle 
correlations. Finally, the Conclusions and scope for future work are presented. 



2 Azimuthal angle correlation ratios in MSYM & QCD 

In this section a review of the observables of interest is presented, following the notation in [17]. 
The kinematical configuration under study is that of Mueller-Navelet jets [29], where two for- 
ward jets with very similar squared transverse momenta p\ and p\ are produced with a large ra- 
pidity separation Y between them, and a relative azimuthal angle d. If a; 1.2 are the fractions of 
longitudinal momentum from the parent hadrons carried by the partons generating the jets then 
Y ~ In (xix%s I P1P2) ■ In the limit of large rapidity separation, terms of the form XY must be 
resummed to all orders, with A being the 't Hooft coupling in MSYM and a s = a s N c /ii in QCD. 
This is equivalent to the Regge limit, where s 3> \jp\p\, which can be treated using the BFKL 
formalism. 

In QCD, for the azimuthal angle correlation ratios to be defined at the end of this section, the 
effect of parton distribution functions is negligible at large Y, therefore the focus will be on partonic 
cross sections. These are written as the convolution of a differential partonic cross section with jet 
vertices $j e t 1>2 : 

&(Pl,P2,Y) = J d 2 q 1 J d 2 q 2 $ jetl {qi , Pi ) d 2^ d 2 q2 ^2 (g2, P2), (1) 
where, for simplicity, the LO jet vertex is used: 

Sjet.faPi) =s agfaPi) = e(q 2 -pi). (2) 

p 2 corresponds to the experimental resolution scale for the transverse momentum of the jet. Con- 
sidering the Green function, /, at NLO, the gluon-gluon differential partonic cross section reads 

da _ ir 2 a 2 f(qi,q 2 ,Y) 



d 2 qid 2 q 2 2 q\q\ 
Using the transformation 



f(qi,q 2 ,Y) = I ^e" Y f( qi ,q 2 ,cj) (4) 
' 2ni 



Dijet Production at Large Rapidity Separation in Af = 4 SYM 



3 



the BFKL integral equation can be expressed in the form 

w/(Qi,Q2,w) = 5 {2) (qi - q 2 ) + J d 2 n K,(q 1 ,K)f(K,q 2 ,u). 
At LO accuracy the solution to this equation is 



f(qi,q 2 ,u) = 



n— — 00 



-00 2tt 2 y/qjqj w - a s xo(|"|, ^) ' 



(5) 



(6) 



where the eigenvalue of the kernel is a function of the logarithmic derivative of the Euler Gamma 
function: 



X o(»,f)=2V(l)-^ 



1 + 71 



+ iv) - ip 



1 + n 



It is convenient to use the following transverse momenta operator notation: 
q\Qi) = 9*1?*); (<?i|l|<?2) = 5 {2} (qi - q 2 ), 



(7) 



(8) 



to introduce the basis of eigenfunctions \n, v) of the LO kernel that appear in eq. (6), satisfying the 
normalization condition (n 1 , v'\n, v) = 8(v — v')8 n ^ n i, as 



1 



{q\n,v) 



nV2 



2 \iv— giwd 



(9) 



At NLO, in the QCD case, this basis is no longer diagonal due to the scale invariance breaking 
stemming from the running of the coupling contributions. The action of the NLO QCD kernel on 
the basis of eq. (9) reads [10] 



(n, v\IC\is' , n) — a 



.-■.MS 



1 



\u I \n'\, - +iv* 



1 



8N C 



d 



l a — % m ~ 21n ' i 
ov av 



s,MS 



A> / d 



(10) 



where a symmetrized version in v and 1/ to express In q 2 as a derivative with respect to either v or 
v' has been used. This is the natural choice in the q\ ~ q\ case, removing terms of the form x'o m 
the kernel, which would be complex for real values of v. All the other functions in eq. (10), with 

v = i ( 5 — 7) , are 



Xi(n, 7 ) = 5xo(«,7) + 2^( 3 ) 



7T 2 C0S(7T7) 



4 sin 2 (^7) (1 - 27) 



7(1 " 7) 
2(3-27X1 + 27)' 



/ N f \ 2 + 37(1-7) \ g f 1+ N f\ 
+ v 1 + iVf J (3 - 2 7 ) (1 + 2 7 ) J " I 1 + iVf J 

n(n, 7) = J \f' (7 + I) + V>" (l - 7 + I) - 20(n, 7) - 20(n, 1 - 7) 
0(n,7)=E 



^'(fc + n + l) -^'(fc + 1) 

^(fe + 1) - + " + 1)" 



u fc+ 7 +f 
+ (-l) fe+1 (/?'(fc + n + 1) + /3'(fc + 1)) + 



(11) 



(12) 



(13) 



fc + 7 + 
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«M7)-*(^)-*(i), *-i («-.> + £). (.4, 

In MS YM the absence of running of the coupling makes the action of the NLO kernel on the LO 
eigenfunctions diagonal and eq. (10) simplifies to 

(n, H^msymK n') = A ( xo (\n'\, I + iiA + A x f SYM (\n'\, I + iA ) 6 n , n ,5(v - i/), (15) 



2 

with A being the non-running analogue of a s , and [10] 2 



xf SYM (n, 7 )- (1 ^ (2)) Xo(|n|,7) + |c(3) + W,7). (16) 

Expressed in terms of the \n, v) basis, the differential cross section in the azimuthal angle 4> = 
i?i — — 7r, with i?j corresponding to each tagged jet, is [17] 



e 



^c„(n d7) 



d^(pf, 2 ,y) = A 2 ~ 

with the coefficients C n (Y) given, in the QCD case, by 

&s(p 2 )Y (xo(|n|,^)+a 8 (p 2 ) (xi(|«|,*)-^ x ( f" U £ 

, , f°° dv e \ \ U+" ) . 

C " (y) " L S TT^ ' <18) 

where, for simplicity, the configuration p 2 ~ p 2 , ~ p 2 has been taken. In the MSYM case, with 
(So = 0, these are 



CJX> = L^- -jt^- - <19) 

Interesting observables can be obtained from these coefficients. C generates the total cross section 
since it is integrated over <fi: 

7r ry 

^( Pl)2! F) = — p=Co(F). (20) 
P1P2 

To study the effect of the higher conformal spins n, which is one of the main points of this work 3 , 
the following ratios of azimuthal angle correlations are introduced: 

<coB(mfl> = (21) 
Cot,*) 

(cos(mfl) _ CrJX) 

m ' n{ } - (cosM)> - C n (Y)- {ZI) 

In the analysis of these quantities different types of coefficients C n will be used. The MSYM co- 
efficients given in eq. (19) will be denoted with C„ . For QCD, when computed with the full 
expression as in eq. (18), C^ LL will be used, keeping C^ L to indicate the use of the LO approxima- 
tion in the exponential. Another useful quantity will be the scale invariant coefficient Cf^ given by 
putting fa = in eq. (18). The 't Hooft limit, defined by N c — > 00 while keeping the product a s N c 
fixed, will be also investigated. 



2 The MSYM calculations in this study have been obtained within the MS scheme, other renormalization approaches give 
very similar results. 

3 The identification of the n's as conformal spins stems from the elastic amplitude, linked to the total cross section via the 
optical theorem, which has an underlying SL(2, C) invariance [30], 
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3 Physical renormalization schemes & BLM procedure 

After the computation of the NLO terms in the BFKL kernel [22, 23], it was noticed that these 
corrections lead to instabilities [31] in collinear regions. All-order resummations of the leading 
collinear contributions helped improve the perturbative convergence [32, 33]. Other proposals based 
on introducing short-range in rapidity correlations were also successful in stabilizing the perturba- 
tive expansion [34, 35]. A complementary approach is that developed in [36, 37] which, instead of 
employing the MS scheme with arbitrary renormalization scale setting, proposed the use of a phys- 
ical renormalization scheme like MOM [25, 26] with optimal scale set by the BLM procedure [28]. 
In this way the NLO corrections have a milder behaviour and the value of the pomeron intercept 
has a very weak dependence on the hard scale of the scattering process, leading to an approximate 
scale invariance. This is a good motivation to compare this approach with MSYM, which enjoys 
four dimensional conformal invariance. 

The transition from the MS to the MOM scheme is equivalent, at LO accuracy, to the redefinition 
of the coupling [25,26] 

«MOM = a MS ( 1 + ^MOM ) , (23) 



where T M om is a function of N c , Nf and of the gauge parameter £: 



rp rriCOnf _i_ rp0 

1 MOM — 1 MOM ' 1 MOM ' 



^mom — 2 I ^ 3^ 

with I = — 2 Jq 1 dxln(x)/[x 2 - x + 1] ~ 2.3439. At NLO this redefinition is equivalent to the 
rescaling 

, _ / TmomA 

/!^/i = /!CXp[ ] • ( 25 ) 

An appropriate choice of renormalization scheme and scale can render the higher order coefficients 
of the perturbative series small. Several prescriptions to do so have been proposed [38-41], with 
possibly the best physically motivated being the BLM one where the coupling redefinition absorbs 
charge renormalization corrections in such a way that the coefficients of the perturbative series are 
identical to those of the corresponding conformally invariant theory with j3 = 0. 

In [36] the BLM approach was applied to the BFKL description of the 7*7* cross section. In 
order to enhance the effect of BLM in gluon dominated processes, it was argued that it is appro- 
priate to use a physical scheme suitable for non-abelian interactions, such as MOM, based on the 
3-gluon vertex [25-27] or the T scheme based on T — > ggg decay. This procedure is not free from 
complications since it has an unnaturally high scale at v = [42]. 

In the present work the BLM setting for dijet production is explored, not only at the pomeron 
intercept level but, very importantly, also investigating its effects in azimuthal angle decorrelations. 
To generalize the BLM scale fixing for arbitrary conformal spins n it is convenient to write the 
BFKL kernel in the form 

1 2 n i . q ms(9 2 )^ ( . . , ,ot m {q 2 )\ 
w ms(<Z > n > v ) =Xo(n,v)— 1 + 7^(71, 1/) . (26) 
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The NLO coefficient r^g is decomposed into /3-dependent and conformal (/3-independent) parts: 



where 



and 



r P {n v) = _?o( xo(n,v) _5\ 



4xo (n,v) 



(^ 2 -4) 



Xo(n, v) - 6C(3) + — sech(7ri/) tanh(?ri/) 
3 2r/ 



71 + 1 



+ W + 0'' 



„ (n + 1 



1 



iz/ — 2(f) I n, — h>f — 2</> ( n, ii/ 



1 



AT, 



1 



A 7 ! / V 4 16(1 + ^ 2 ) 



S n ,o ~ 1 + 



1 



y V8 32(l + ^ 2 ) 



(27) 



(28) 



(29) 



The NLO BFKL intercept in the MOM-scheme, evaluated at the optimal BLM scale can be repre- 
sented as follows: 



n ( a 2M0M )N I a ffr 2MOM ) 



(30) 



where r M 0M(", f) = r Ms( n ; ^) + ^mom- The BLM redefinition of the coupling involves the (3- 
dependent part of the corrections implying that 



S 2 BTM>^)=<7 2 exp 



MOM 



4r MOM( n ^) 



2 1 i 1+4/ 



(31) 



with r^ OM (n, i/) = r^(n, v) + T MO m- 



4 Comparing MSYM with QCD results 

As a first result, in Fig. 1 the intercept of eq. (30) for n = computed in the MOM scheme with 
BLM scale is confronted with the MSYM intercept and the results at LO and NLO with no BLM 
scale fixing. The coupling for MSYM is chosen in a range between A = <5 s (q 2 /4) (MSYM_) 
and a s (4q 2 ) (MSYM + ), which corresponds to the light yellow band. This plot is in agreement 
with the calculation performed for the 7*7* total cross section in [36] since the scale invariance, 
with respect to the photon virtualities in that case and the jet transverse momentum now, of the 
asymptotic intercept in the MOM scheme is manifest. This intercept for the conformal invariant 
MSYM theory at NLO is very close to the LO one, already hinting towards a better perturbative 
convergence than in QCD. For completeness, results for the QCD calculation taking the 't Hooft 
limit N c — > 00 are also shown. It is important to note that the MOM-BLM is the closest to MSYM 
of all renormalization schemes in QCD, in explicit agreement with the idea of it being collecting the 
conformal contributions to the observable. 

In FlG. 2 a typical scale p 2 = 15 GeV 2 is taken to show, for QCD, the asymptotic intercepts for 
each of the azimuthal angle Fourier components n in the MOM-BLM scheme, which can be read 
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from the values of the curves at v — 0. The dominant component at large relative rapidity separation 
of the dijets is n = with all the other intercepts being negative. A similar behaviour is found for 
other schemes and in the MSYM case. In the following it will be investigated how these intercepts 
affect physical observables sensitive to the azimuthal angle dependence. 

0.6 
0.4 

0.2 


-0.2 

o -0-4 
UJ( P 2 ,0) 



-0.8 
-1 
-1.2 

10 20 30 40 50 60 70 80 90 100 

9 

P 

Figure 1: Comparison of the evolution of the intercept with the value of the jet resolution scale p 2 
in LO and NLO for different renormalization schemes in QCD and the MSYM case. The curve 
W'tHooft(p 2 ) corresponds to taking the 't Hooft limit in the QCD NLO case. 

As it was discussed in previous sections, the azimuthal angle dependence of the gluon Green 
function can be studied by expanding it in Fourier components, each of them characterized by a 
discrete parameter n, which in elastic scattering has the interpretation of a conformal spin associated 
to an underlying SL(2, C) symmetry. In this way, by studying observables sensitive to the azimuthal 
angle, the Moebius invariance present in gauge theories at high energies is under scrutiny. 

As eq. (20) shows, the n = coefficient is directly related to the normalized cross section: 
a(Y)/a(0) — Co(Y) /Co(0). The rise of Co with the rapidity separation between the two tagged jets 
is shown in Fig. 3. The first feature to be noticed is the faster growth of the MSYM cross section 
for the wide range of values in the 't Hooft coupling, A, considered in the yellow band. This shows 
that the NLO real emission kernel in the MSYM theory dominates over the virtual contributions in 
a much stronger fashion than in the QCD case, for any renormalization scheme. This also indicates 
that the net effect of introducing the extra fields in the supersymmetric multiplet increases the minijet 
multiplicity in the final state. In future works it will be worth looking into these details of the final 
state using event generator Monte Carlo techniques [43]. 

It is noteworthy that for small rapidity separations the QCD cross section in the MS-scheme is 
lower than in the MOM-BLM scheme, with the latter being closer to the MSYM result. This is 
to be expected from a renormalization set-up which resums the higher order dominant conformal 
contributions. However, at Y ~ 10 a departure from this behaviour is observed, with the calculation 
in the MS-scheme being now closer to the MSYM one. This again hints at the instability of the 
n = conformal spin component, which is known to be very sensitive to collinear radiation, not 
fully included in the BFKL kernel. It is then natural to predict a similar crossing of behaviour at 
some Y for any quantity sensitive to the n = eigenvalue of the kernel. This is indeed what is found 
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Figure 2: Eigenvalues of the QCD BFKL kernel at NLO in MOM-BLM scheme, for different con- 
formal spins n. 
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Figure 3: Growth with dijet rapidity separation of the cross section in MSYM and QCD for different 
renormalization schemes. 
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Figure 4: Evolution of the average of cos (j> with jet rapidity separation in MS YM and QCD for 
different renormalization schemes. 
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Figure 5: Evolution of the average of cos(20) with jet rapidity separation in MSYM and QCD for 
different renormalization schemes. 
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when the average of cos {nip) is calculated as in eq. (21). A couple of examples with a cross-over 
of lines are plotted in Figs. 4 and 5 for the n = 1, 2 cases. In these figures it is also interesting 
to note that the MSYM dijets are less correlated in azimuthal angle in the MSYM theory than in 
QCD, which corresponds to another manifestation of the higher multiplicity of partem radiation in 
the supersymmetric case. 

1 



K 2 ,i 01 



0.01 

2 4 6 8 10 12 14 

Y 

Figure 6: Evolution of 72-2,1 = ^f0j^r with jet rapidity separation in MSYM and QCD for different 
renormalization schemes. 

The main conclusion of this analysis is that to define observables only sensitive to conformal 
dynamics it is needed to remove the n = contribution. One way of doing this is to use the ratios 
of azimuthal angle averages H m ^ n (Y) in eq. (22). Two of these ratios (72-2, i> ^3,2) are calculated 
in Figs. 6 and 7. It is important to note that the MSYM ratios are very close to those calculated 
in QCD, indicating that these observables capture the bulk of the conformal dynamics in QCD. 
Moreover, the renormalization scheme which gives ratios closest to those in MSYM is the MOM- 
BLM scheme, independently of the separation in rapidity between the two tagged jets. Having 
removed the n = dependence, the cross-over of lines does not take place anymore. 

5 Conclusions 

In this work dijet production has been investigated for configurations where the two tagged jets 
are largely separated in rapidity. This has been done both in QCD and MSYM, with particular 
emphasis in investigating how different renormalization procedures in QCD best reproduce the 
conformal dynamics of MSYM. With this purpose, ratios of azimuthal angle decorrelations, 
which have an excellent perturbative convergence, have been shown to capture the bulk of the 
conformal contributions. For these ratios the results calculated in QCD with the MOM-BLM 
scheme are very similar to those obtained in MSYM. It has also been shown that the two 
tagged jets are less correlated in azimuthal angle in MSYM than in QCD, indicating that in the 
supersymmetric case there is a larger multiplicity in the final state. For future work it will be inter- 
esting to investigate other multijet configurations and different observables at a more exclusive level. 
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